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We investigate Laughlin's fractional quantum Hall effect wave function in the cylinder geometry of Laugh- 
lin's integer quantum Hall effect argument, at filling factor 1 /3. We show that the plasma analogy leads to 
a periodic density, and that the wave function admits a representation as a "quantum polymer", reminiscent 
of the quantum dimer model by Rokhsar and Kivelson. We explain how the representation can be exploited 
to compute the normalization and one-particle density in the limit of infinitely many particles. 



1 Introduction 

Laughlin's wave function [1] is widely accepted as a good description of fractional Hall effect ground 
states at simple filling fractions. It was initially proposed for a disk geometry, but later adapted to various 
geometries [2-A]. In this article, we study Laughlin's function on a cylinder [2]. The choice of geometry 
is motivated by its use in Laughlin's argument [5] for the integer quantum Hall effect (IQHE) and the 
suggestion that a ground state degeneracy is required in order to reconcile Laughlin's IQHE argument with 
fractional charge transport [6]. 

Ground state degeneracy, in turn, is closely related to broken translational symmetries and multiple 
Aharonov-Bohm periods (see, e.g., [7, 8], and Sec. 2). In an earlier paper [9], we have shown that Laugh- 
lin's state on thin cylinders, at filling factor 1 /3, is periodic, with period three times that of the filled Landau 
level. In this article, we present considerations complementary to [9]. 

A first aspect pertains to the plasma analogy, which relates the modulus squared of Laughlin's wave 
function to the Boltzmann weight of a classical two-dimensional one-component plasma, or jellium. It was 
used in [1] to justify that Laughlin's wave function in the disk, at filling fraction 1/3, has a homogeneous 
density. In Sec. 3, we explain that the same analogy, as soon as it is applied to the cylinder, leads to a 
periodic density. This uses results on jellium on cylinders ("jellium tubes") [10]. 

A second aspect focuses on a detailed analysis of the structure of Laughlin's wave function derived 
from combinatorial properties of powers of Vandermonde determinants [11,12]. We show that Laughlin's 
wave function admits a representation that we call a "quantum polymer" because of its resemblance with 
the quantum dimer model [13] (Sec. 4). It turns out that this representation is very useful to access the 
normalization and the one-particle density: In Sees. 5 and 6, we review the results on normalization and 
the density's periodicity obtained in [9] and propose a computational scheme, based on the polymer repre- 
sentation, for the normalization and the density in the limit of infinitely many particles. Roughly, we obtain 
an expansion in powers of cxp(— j 2 ), where 7 is the ratio of the magnetic length and the cylinder radius. 
The zeroth order in the expansion is the Tao-Thouless state [14], which is 3-periodic and corresponds to the 
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limit of infinitely thin cylinders. This is closely related to results by Rezayi and Haldane [15]. However, 
periodicity does not stop at order zero. In fact, we obtain a sequence of approximate densities, each of 
which is periodic with smallest period 3. 

Hence, in the FQHE context, the quantum polymer representation is useful in two ways. First, it shows 
that a 37-periodicity is indeed built-in into Laughlin's wave function on a cylinder, as is expected [6] from 
Laughlin's IQHE argument. Second, the quantum polymer formalism gives a practical way to access the 
one-particle density of Laughlin's function. Our method can actually be extended to handle correlation 
functions. This is interesting in view of the widespread use of Laughlin's function in theoretical investiga- 
tions of FQHE ground states. 

2 Symmetry breaking, ground state degeneracy and multiple Aharonov- 
Bohm periods 

Before we turn to the analysis of Laughlin's wave function, let us briefly recall the relationship between 
symmetry breaking, ground state degeneracy and multiple Aharonov-Bohm periods in the cylinder geom- 
etry of Laughlin's IQHE article [5]. 




(3JV - 3)7 



Fig. 1 Cylinder geometry and lowest Landau level basis functions. The basis functions are Gaussians 
with centers spaced apart by 7. Laughlin's state occupies a cylinder of length roughly (37V — 3)7. 

Consider electrons moving on a cylinder of radius R in a perpendicular magnetic field of strength B (see 
Fig. 1). The cylinder is threaded by an additional flux ip, which does not change the magnetic field B but 
changes the vector potential. Each particle position is described by its coordinate x along the cylinder axis, 
and by an angular coordinate y 6 [0, 2tvR[. It is convenient to choose units such that Planck's constant 
H, the electron mass m e , the elementary charge e and the field magnitude B all take the value 1. In those 
units, the magnetic length i = (h/eB) 1 I 2 equals 1 and the flux quantum is <po = h/e = 2n. A crucial role 
is played by the dimensionless parameter 7 = l/R which in our units becomes the inverse of the radius, 
7 = 1/R. In the absence of a background potential, the one-particle Hamiltonian for the infinite cylinder 
is 

H^) = \[-dl + {-\d y -x+^-f] (1) 

with periodic boundary conditions in the y-direction. 

The family of Hamiltonians H (cp) has two intimately related features, gauge periodicity and a discrete 
translational invariance, see Fig. 2. Both aspects become explicit in the lowest Landau level, i.e., the 
ground state of the Hamiltonian (1). It has a complete orthonormal set 

ipk{z;<p)= , = =cxp(ifc7y)exp{-^[x- (fc + y/27r)7] 2 }, z = x + iy, k G Z. 

Notice that \ipk(z;(p)\ 2 , k e Z, are Gaussians of x with centers spaced apart by 7. Consequently, the den- 
sity in the filled Landau level is a sum of equally weighted Gaussians. Its Fourier series can be computed 
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exp(— iwypx] 

usual transl. 
n real 



H(tp + 27Ttt) 



exp(in7y) 


gauge transf. 
n integer 



t(nyex) — expfinyy) cxp(— ip^rry) 
magnetic transl. 
n integer 



Fig. 2 Gauge periodicity and discrete trans- 
lational invariance for the Landau Hamilto- 
nian (1) on a cylinder: a usual translation by 
717 changes the flux from ip to tp + 2nn (full ar- 
row). When n is an integer (dotted arrows), the 
change in flux can be gauged away, leading to 
the invariance of the Hamiltonian with respect 
to the magnetic translation. 



using Poisson's summation formula, 

CO 

PfiiicdLLL(z; v?) = ^2 Hk(z;<p)\ 2 

k— — oo 

This formula shows that the filled Landau level density is always periodic with minimal period 7, although 
the amplitude of oscillations becomes very small when 7 is small. 1 Changing if amounts to an overall 
shift of the density. After a change of <p by 2w, the density returns to itself: the translational period 7 is 
associated with the simple Aharonov-Bohm period Aip = 2tt. 

Now suppose that a ground state has translational period 37 instead of 7. Then it takes 3 flux units before 
the state returns to itself. Because of the system's gauge periodicity, the states obtained after addition of one 
or two flux units are ground states as well. Hence, 37-periodicity leads not only to a threefold Aharonov 
Bohm period but also to threefold degeneracy of ground states. 



= ^l 1 + 2£exp(— — )cos[fc( — - ¥>)]}. 



3 Plasma analogy 

Laughlin invoked the plasma analogy to show that the density of his wave function, in a disk ("Laughlin 
droplet"), is constant. In this section, we show that the same analogy, when it is invoked in a cylinder ge- 
ometry, leads to a one-particle density that is periodic with period 37. This is based on earlier observations 
that jellium tubes should have a periodic density [10, 17]. For plasma parameter r = 2, this was already 
shown in [18]. 

Let us start by a brief description of the system we are interested in. We consider N particles of 
charge q moving on a strip [—L/2, L/2] x [0, 27ri?] and interacting with a neutralizing background of 
uniform charge density — nq. We work with complex coordinates z = x + iy and choose a semi-periodic 
logarithmic "Coulomb" potential V(z) := — In |2 sinh[z/ (2i?)]|. The potential energy is 

N pL/2 p2-!rR 

<5>{z 1: ...,z N ) = q 2 V( Zj - z k )-nq 2 ^2 / V(zj - z)dxdy 

l<j<k<N j=l J~ L / 2 J® 

+ \n 2 q 2 [ [ V(z- z')dxdydx'dy'. 

2 J[-L/2,L/2] 2 J[0,2irR] 2 

For the special choice that the plasma parameter T := (3q 2 equals twice the exponent in Laughlin's func- 
tion, here, r = 2 • 3, the Boltzmann weight cxp(— /?$) becomes proportional to Laughlin's function. More 
precisely, let us choose units such that the background density is n = (3 ■ 27T)" 1 , and suppose that V = 6. 
Then 

3(7Y-1) 



cxp(-/?$(zi, ...,z N )) oc $ N {{Zj 2 7 J') 



(2) 



Here, ^jv is Laughlin's cylinder function with magnetic length 1=1, see Eq. (3). Therefore the plasma 
system and Laughlin's function will have the same one-particle density. 



The amplitude is, roughly, cxp(— 7r 2 /7 2 ) = exp(— tt 2 R 2 ). A similar small 7 behavior for the amplitude of oscillations of 
the 37-periodic density at filling factor 1/3 has been found numerically in [16]. 
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For a plasma on a plane, at coupling T = 6, we expect a homogeneous density. On the other hand, 
a cylinder is more like a quasi one-dimensional structure, especially when we consider long cylinders 
[17]. Notice also that the y average of the potential V(z) equals —\x\/(2R), i.e., a one-dimensional 
Coulomb potential. As was observed in [10], the comparison with one-dimensional jellium leads us to 
the expectation that the density is periodic in the direction of the cylinder axis with period L/N. This 
is because one-dimensional jellium always displays symmetry breaking [19,20]. The interpretation is 
"Wigner crystallization", the electrons tend to minimize their repulsive interaction by adopting positions 
on a lattice with spacing L/N. 

Now, the crucial point is that the period L/N expected for jellium on a cylinder translates into 37- 
periodicity for Laughlin's function on a cylinder. Indeed, the cylinder's length is 



2nRn 2TTJ- 1 (3 ■ 271-)- 1 ' 

Thus because of the plasma analogy, we expect that the density in Laughlin's function is periodic with 
period 37, i.e., three times the period of the filled Landau level. The periodicity is confirmed by numerical 
results on jellium tubes [10], from which it appears that the amplitude of oscillations is small when the 
cylinder radius gets large. This fits nicely into the picture that the cylinder interpolates between a line and 
a plane. The important point is that the cylinder retains the periodicity of the one-dimensional system. 



4 Laughlin's wave function as a quantum polymer 

In the previous section, we have used the plasma analogy to argue that a 37-periodicity in the axial direction 
is very natural. In this section, we take a different approach and investigate directly the structure of the wave 
function. This leads us to a recursive pattern in the sequence of Laughlin functions (^n) n^n, formalized 
with a "product rule", and to the representation of Laughlin's function as a "quantum polymer". 

We use the units and notation of Sec. 2, with the flux parameter ip set to 0. Laughlin's function adapted 
to the cylinder geometry [2] at filling factor 1/3 is 

N 

#jv(zi, — ,zn) = kn (exp(~fZk) - exp(~/Zj)) exp(- y^/2). (3) 

l<j<fc<JV j=l 

It will be convenient to fix the multiplicative constant kn as 

11 9 N ^ 

KN = TW a^FwW^ 2 exp( ~2 72 g n (4) 



4. 1 Product rule 

The wave function is a polynomial of exp(7Zi), exp^z^v) times a Gaussian weight. Expanding the 
polynomial, we find an expression of "f^ in terms of lowest Landau level basis functions tpk(z) — ipk(z; 0) 
given in Sec. 2. For up to three particles, this gives 

*i = Vk>, 

*2 = "00 A tp 3 - 3e" 27 Vi A -02 

*3 = i' A iJj 3 Alps — 3e~ 27 V'o A ip4 A ips - Se^^i'i A ip 2 A ?/>6 
- Qe^tpi A -03 A -05 + Ibe' 8 ^^ A ^3 A ip A . 
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Notice that the coefficient —3 exp(— 27 2 ) in the expansion of ^> 2 reappears in the expansion of ^3. Simi- 
larly, the coefficients in the expansion of \& 3 will show up in the expansion of ^4. The aim of this section 
is to explain and formalize this phenomenon. 
In general, one can write down an expansion 

= ^2 a N {n%i, ...,m N )ip mi A ... A tp mN . (6) 

0<m 1 <...<m N <3N-3 

As can be seen in Eqs. (5), the expansion coefficients awim) are always made up of two parts, an integer 
and a power of exp(— r y 2 ). The integer comes from the expansion of the third power of a Vandermonde 
determinant. The power of cxp(— j 2 ) accounts for the multiplicative constant k n and for the normalization 

of cxp(jkz) exp(— x 2 /2). More precisely, 

1 9 N 

a N (mi, ...,m N ) = b N (mi, m N ) exp(-7 2 ^(m 2 - 9(j - l) 2 )) 

with suitable integers (m) . The coefficients (m) are specific to the cylinder geometry and have been 
analyzed by Rezayi and Haldane [15], while the integers bpf(m) play a role as well for the usual "Laughlin 
droplet" geometry and have inspired a certain amount of combinatorial research [11, 12, 21, 22]. 

It is practical to associate with each finite sequence m a sequence of l's and O's characterizing occu- 
pancy or non-occupancy of lattice sites { — 1, 0, 1, 37V — 2}. Each site k stands for the orbital ijjk- The 
sites —1 and 3iV — 2 will never be occupied, but it is advantageous to include them nevertheless. For 
N = 1, there is only one configuration, 010, corresponding to mi = 0. For N = 2, the two wedge 
products give rise to 

(0,3) (010 I 010), (1,2) <-> (001100). 

Observe that (0, 3) corresponds to a concatenation of two identical one-particle blocks 010, as highlighted 
by the vertical bar. Similarly, for N = 3, one- and two-particle blocks show up: 

(0,3,6) ^ (010 I 010 I 010), (1,2,6)- (001100 , 010), 
(2,3,4) ^ (000111000), (0,4,5) (010 | 001100), 

etc. Now, we assign to each block the amplitude of the corresponding m. Thus 

ai(010) = ai ( mi = 0) = 1, a 2 (001100) = a 2 (l, 2) = -3e" 272 . 
The crucial observation is the following product rule: 

Concatenation of blocks results in the multiplication of amplitudes. 
For example, 

a 3 (010 I 010 I 010) = a 3 (0,3,6) = 1 = a!(010)ai(010)ai(010), 

(8) 

a 3 (010 I 001100) = a 3 (0,4,5) = -3e~ 27 = ai (010)a 2 (001100). 

A rigorous proof of the product rule for aAr(m) is given in [9], the product rule for the integer part &Ar(m) 
has been shown earlier [12]. 

There is a simple characterization of the block decomposition of a sequence of occupation numbers 
corresponding to a vector m in terms of the differences 



k k 

^:=^m,-^3(i-l), ke{l,...,N}. 

3=1 i=i 



(9) 
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N(X) X u x a N (X) = \\ux\\' 2 

1 {-1,0,1} 1 

{2,3,4} V>3 1 

{5,6,7} 1 

~~ 2 {-1,...,4} -3e- 2 T>iA^ 2 9e-^ 

{2, ...,7} -3e- 2 ^ 2 tj; 4 Atjj 5 9e~ 4 ^ 2 

~~ 3 {-1,...,7} -6e- 5 T>i A tl> z A fa 36e- 10 ^ + 225 e - 16 ^ 



+15e~ 87 ^2 A^A^ 

Table 1 Polymer functions ujf contributing to the representation (10) of Laughlin's wave functions for 
N = 3 particles. 

When m has a non-vanishing amplitude ajy(m), the i/j-'s are non-negative and = [12]. The block 
decomposition is determined by the additional zeros of (i>k)'- each bar after 3fc sites corresponds to a 
vanishing difference, i.e., Vk = 0. This is the characterization used in [9, 12]. 

A different characterization is in terms of squeezing operations [15]. Rezayi and Haldane [15] observed 
that every m with non vanishing amplitude apf(m) can be obtained from (0, 3, 6, 3(A — 1)) by a series 
of elementary operations (oi, a/v) i— > (a%, dj + 1, dj — 1, ajv), i < j, or "squeezings". If for 
some k we only consider sqeezings that respect a boundary at k, i.e., never squeeze pairs with i < k 
and k < j, then the resulting m will split into two blocks of k and N — k particles. 

It is important to notice that the product rule contains, implicitly, a periodicity statement: the amplitude 
assigned to a block does not depend on its position in the sequence of occupation numbers. In Eq. (8), 010 
has amplitude 1 regardless whether it refers to a particle sitting in the site 0, 3, or 6. This periodicity is at 
the origin of 37-periodicity in Laughlin's state. 

Remark 4.1 (Reversal invariance) In addition to the product rule, there is another property of the 
coefficients 6jv(m), shown in [11,12], that transfers to the amplitudes aN(m): replacing m with the 
reversed sequence (3./V — 3 — mjv-j)j leaves the amplitude unchanged. This means that Laughlin's cylinder 
function is unchanged by overall rotations by 180° around the middle of the cylinder, (3./V — 3)7/2. 

4.2 Quantum polymer 

Now we translate the product rule for the amplitudes ajy (m) into a representation formula for the function 
\&jv. As a first step, notice that the way a configuration splits into blocks defines a partition of the discrete 
volume Ajy = { — 1, 0, 3N — 2} into discrete intervals, or "rods", Xj. For example, 

m = (0, 3) = (010 I 010) ^ A 2 = {-1, 0, 1} U {2, 3, 4} = X 1 U X 2 . 

In the expansion (6), we can group together m's that belong to the same partition. Using the product rule, 
we obtain a representation of ^jy as a sum over partitions. For N = 3 particles, this gives 

*3 = ^{-1,0,1} A 1i{ 2 ,3,4} A "{5,6,7} 

+ "{-1,0,1} A W{2,...,7} + "{-1,...,4} A lt{ 5 ,6,7} + U{-1,...,7} 

with the functions ux from Table 1. More generally, with each rod X of cardinality (or "length") a 
multiple of 3, we associate a fermionic wave function ux of |-X"|/3 complex variables, ux describes a 
cloud of electrons of density 1/3, localized in the discrete volume X: 

u x G A N ( x hpim{i>k I k G X}, N(X) = \X\/3. 
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Shifting a rod amounts to a magnetic translation of the function: 

U3 j+ x = iO'3 7 e x ) 0Af(x V, 3 e Z. (11) 

This translational covariance is a consequence of the periodicity in the amplitudes, mentioned in the pre- 
vious subsection. With the functions ux, Eq. (10) is generalized as follows: 

Laughlin's wave function ^ n at filling factor 1 /3 is a sum over ordered partitions of the discrete 
volume {— 1, 0, 3-/V — 2} into rods Xj with lengths that are multiples o/3: 

*at= u Xi A ... A ux D ■ (12) 

(X U ...,X D ) 

A partition is ordered when the rods Xj are labelled from left to right. In Eq. (12), different partitions give 
actually rise to orthogonal contributions. 

The representation (12) is reminiscent of the quantum dimer (QD) model proposed by Rokhsar and 
Kivelson [13] as an idealization to resonating valence bond states. QD wave functions are given as sums 
over partitions of a lattice into dimers (i.e., subsets of two neighbor elements). Each dimer represents two 
spin 1/2 particles that form a singlet. Motivated by the quantum dimer model, we will say that Laughlin's 
function can be expressed as a (discrete) quantum polymer. The associated polymer system consists of 
3n-mers and has a built-in periodicity, see Eq. (11). 

In the next sections, we will exploit the quantum polymer representation to obtain information on the 
normalization and the density of Laughlin's function. 



5 Normalization 

For computational purpose it is useful to have information on the normalization constant Cn = H^atII 2 
of the wave function. In [12], the plasma analogy is used to derive results on the asymptotics of the 
normalization in a disk geometry. Here, a similar approach shows that 

lim ^-logCjv = -logr 

N — >oc TV 

for some r > 0. The limit is closely related to the free energy of the associated classical plasma, see Sec. 3. 
The quantum polymer representation allows us to give more refined results, namely the convergence 

lim r N C N = q>0. 

JV->oc 

The quantities r and q are determined by the polymer functions (ux ) in a simple way. Because of the 
translational covariance (11), the norms \ \ux\\ depend on N(X) only and we can define non-negative 
numbers (a n ) n >i by 

otN(x) ■= \\ux\\ 2 , 

see also Table 1 . The polymer representation leads to a formula for the normalization, 

C N = ||*jv|| 2 = ^ IKyJI 2 • ■■■ • ||wx D || 2 = 2J a ni -...-a 7lD . (13) 

(X!,...,X D ) D, ni ,...,n D :£N 

■rn + ...+n D = N 

It follows that (Cn) satisfies a recurrence known in stochastics as a discrete renewal equation [23], with 
coefficients (a n ) (see Eq. (15)). As a consequence, one of the following cases necessarily holds: 

!• X)^Li r ™ a ™ = l,fJ,:= nr n a n < oo, q = pT x > 
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2 - X^=i rna n = 1, A* := D^Li nr n a„ = oo, q = 

3- En=i < 1. 1 = 0- 

In the first two cases, r is the unique solution of the equation t n a n = 1. In the third case, this equation 
has no solution and J2 n t n OL n diverges as soon as t > r. 

A probabilistic interpretation can be found in [23]. We give an interpretation in terms of the associated 
polymer system. Eq. (13) says Cjy is a polymer partition function with activities (a n ), see [24]. In this 
context, the rescaling of activity in the way a n i— > r n a n is very natural, since it leaves polymer correlation 
functions unchanged - much in the same way as quantum-mechanical expectation values (iff, A^f) /| |\&|| 2 
are unaffected by the multiplication of "J with an overall factor, $ — > c^> . The previous cases may be 
described as follows: 

1. The activities (a„) rescale to a probablity distribution on N. The expected polymer length is finite. 

2. The activities (a n ) rescale to a probablity distribution on N with infinite expectation value. 

3. The activities cannot be rescaled to a probability distribution. 

The last two cases are slightly pathological, and we would like to exclude them. This can actually be done, 
for large 7 (thin cylinders), by a perturbative argument. What happens for small 7 remains open. 

Rezayi and Haldane [15] have observed that for a fixed, finite number of particles, the limit 7 — > 00 
of Laughlin's wave function gives the Tao-Thouless state; in block notation, this is a pure trimer state 
(010 I 010 I ... I 010). In the polymer picture, one can show that the activities satisfy 

ai = l, on = 0(e- 4( - n -^ 2 ) (7 ->(»), (14) 

see also Table 1. Hence, in the limit of infinitely thin cylinders (7 — > 00), all activities except ati vanish. 
In this limit, Cn = l,r= 1, and we are in the first case, i.e., q = lim r^CV > 0. One can show that this 
survives for sufficiently large 7 (see [9] for details). 

A computational scheme. Closely related to Eq. (14) is the observation that long polymers have a small 
activity. This suggests approximating Laughlin's state by a quantum polymer with polymers that cannot 
become longer than 3m, for some cutoff parameter m E N. The normalization of such a cutoff polymer 
system satisfies a recurrence relation of finite order, and the normalization can be computed explicitly. 

We are thus let to the following procedure for computing approximations r m , q m to the quantity r, q 
that appears in the asymptotics of (Cjv)- 

• Compute the normalization of Laughlin's function for up to m particles, Ci, C m . 

• Compute the activities ot\, a m of polymers of length smaller or equal to 3m, by using the recur- 
rence relation 

C n =a 1 C n - 1 + ... + a n - 1 C 1 +a n (n > 2). (15) 

• Solve the equation ait + o^i 2 + a^t 3 + ... + a m t m = 1, t > 0. The unique solution r m is an 
approximation to r. An approximate value to q = liniAr^oo r N Cn is given by q m := J^T n0L rJ' r !n- 

For each fixed 7, the approximations r m , q m converge to r, q when the cutoff parameter m goes to infinity. 
For sufficiently large 7, we can say something about the speed of convergence: 

r m - r = 0( e - 4 (™" 1 )T 2 ) ( m -> 00) 

and similarly for q m . 
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Roughly, this procedure is a way of computing the expansion of r in powers of cxp(— j 2 ). The inter- 
esting point is that each approximate value r m only requires knowledge of Laughlin's wave function for 
a finite number of particles. Notice also that in this scheme the Tao-Thouless state, with normalization 1, 
corresponds to the first approximation, m = 1 (only trimers allowed). 

In the next section, we will see that something similar can be done for the one-particle density and the 
occupation numbers. 

6 One-particle density 

Now we explain how the polymer representation and the results on the normalization lead to symmetry 
breaking, and how the periodic density can be computed. 

Recall that the one-particle density pn(z) is obtained from |^jv(zi, zm)\ 2 by integrating out all 
variables but one, and normalizing in such a way that J piy(z)dxdy = N. As usual, let [-^/2J denote the 
largest integer below N/2. In order to eliminate boundary effects, we take the limit of long cylinders (fixed 
radius), and obtain a periodic density: 

For sufficiently large 7, the one-particle density Pn(z) of Laughlin's wave function (3) has a 
limit p(z) when N — > 00 for fixed radius R = I/7: 

p(z) = lim p N (z-3[N/2\j). 

N—>oc 

The limiting density p(z) is periodic in the direction along the cylinder axis. The smallest period 
is 37. 

The shift 3\N/2\^ ensures that we look at the density near the middle of the cylinder, far away from 
boundaries. - A formal proof can be found in [9]. Here, we shall focus on the resulting representation of 
the periodic density p{z). 

Let fik be the number operator for the lowest Landau level basis function tpk(z). In the limit N — > 00, 
the density p(z) is a sum of Gaussians weighted by the occupation numbers (fik)'- 

00 00 
p(z)= £ (h k )\Mz)\ 2 * £ (n fe )cxp[-(z-fc 7 ) 2 ]. (i 6) 

k— — oo k— — oo 

A similar formula holds for pn{z) (finite N); in this case the summation goes from to 3iV — 3. 

When 7 is large enough so that the limit q of r N CV is strictly positive, there is a nice formula for the 
infinite cylinder occupation numbers (fik)' 

("k) = ^2p F ( x )(nk)x, (fik)x = (ux,fikUx)/\\u x \\ 2 - (17) 
x 

The sum is over polymers 

X = {3j-l,3j,...,3(j + n-l) + l}, jeZ, neN, 
and the "polymer correlation" 

p p (X) ~qr N ^a N{x) (18) 

is interpreted as the probability of finding the rod X. The number (fik)x may be thought of as the proba- 
bility of finding a particle in the lattice site k given that k is covered by the polymer X. 

Eqs. (17) and (18) are proved in [9]. Here, let us only mention that p F is uniquely determined by the 
following criteria: 
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• It is 3-periodic: p p (X + 3) = p F (X). This reflects the covariance (1 1). 

• The probability for finding a polymer X is proportional to the probability r ff ^' aiv(x) (recall from 
Sec. 5 that (r n a n ) defines a probability measure on N). 

• The polymer correlation p F (X) satisfies the sum rule J^xso P F (X) = 1. The sum rules expresses 
that each lattice site (in particular, 0) is covered by exactly one polymer [24]. 

Through Eq. (17), the occupation numbers and the one-particle density inherit the 3-periodicity (11) 
from the polymer system: 

(nfe +3 ) = (fik), p{z - 37) = p(z). 

Thus 3 (resp. 37) is one period; but is it the smallest period? 

The answer is yes. Again, the proof goes with a perturbative argument. Because of Eq. (14), only trimer 
blocks "010" survive in the limit 7 — > 00. As a consequence, [9] 

K>=('r o( r p «:7 2)) ' " £3Z ' <>■» 

I C(cxp(— 7") J , else. 

This generalizes a similar statement on the occupation numbers for a fixed, finite number of particles [15]. 
It follows from Eq. (19) that for large 7, the smallest period of p(z) is 37. Numerical results [10] show that 
the period is 37 for all values of 7 > 0, although the amplitude of oscillations becomes very small when 7 
is small. 

Remark 6.1 (Average density, reversal invariance) With our formulas one can check that the one- 
particle density oscillates around the average value (3 • 27r)~ 1 = (3 • 2ir£ 2 )~ 1 , which is the correct value 
for filling factor 1/3. - The reversal invariance (see Rem. 4.1) leads to a a symmetry with respect to the 
origin: we have (n_fc) = (rife), p{—z) = p(z) (see also [10]). 



An approximation scheme for occupation numbers and the one-particle density. In Sec. 5, we 
have explained how we can compute approximations to the quantities involved in the asymptotics of the 
normalization constants by cutting off long polymers. Something similar can be done for the occupation 
numbers (rife). Notice that, because of Eq. (16), computing occupation numbers and the one-particle 
density is essentially the same. 

Let m € N be a cutoff parameter as in Sec. 5. We consider the quantum polymer system obtained 
from Laughlin's state by discarding all polymers with length strictly greater than 3m. In the limit of long 
cylinders, the one-particle density of this truncated system is periodic, regardless of the size of 7. The 
occupation numbers are obtained by modifying Eqs. (17) and (18) as follows: In Eq. (17), the summation 
is restricted to rods X with length \X\ < 3m, and in Eq. (18), r and q are replaced with their truncated 
values r m and q m as defined at the end of Sec. 5. 

In this way, for each m and every fixed 7, we obtain approximate occupation numbers (hk)^" 1 '- In 
order to compute the m'th approximation, we need only know Laughlin's wave function for up to m 
particles. Again, the first approximation (m = 1) is just the Tao-Thouless state with occupation numbers 
1 at multiples of 3 and otherwise. Approximations for larger m look more complicated, but for each m, 
the approximate occupation numbers are 3-periodic. 

It remains to see if this approximation scheme converges. For the normalization, we know that r m — > r, 
qm — ► q for each fixed cylinder radius. The answer for the occupation numbers depends on the validity 
of Eq. (17): so far we only know that Eq. (17) holds when 7 is large enough. In this case, the approx- 
imate occupation numbers do converge to the occupation numbers in Laughlin's state, and the speed of 
convergence is 

(n k ) [m) - (rife) = 0(e- A( - m -^ 2 ) (m -> 00), 
which complements Eq. (19). 
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